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Abstract 



We study the coherent dynamics of globally coupled maps showing macro- 
scopic chaos. With this term we indicate the hydro dynamical-like irregular 
behaviour of some global observables, with typical times much longer than 
the times related to the evolution of the single (or microscopic) elements of 
the system. The usual Lyapunov exponent is not able to capture the essen- 
tial features of this macroscopic phenomenon. Using the recently introduced 
notion of finite size Lyapunov exponent, we characterize, in a consistent way, 
these macroscopic behaviours. Basically, at small values of the perturbation 
we recover the usual (microscopic) Lyapunov exponent, while at larger val- 
ues a sort of macroscopic Lyapunov exponent emerges, which can be much 
smaller than the former. A quantitative characterization of the chaotic motion 
at hydro dynamical level is then possible, even in the absence of the explicit 
equations for the time evolution of the macroscopic observables. 

PACS: 05.45.+b 

KEYWORDS: Macroscopically coherent dynamics, Finite size Lyapunov ex- 
ponent, Globally coupled maps 
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I. INTRODUCTION 



In the last years the emergence of non trivial collective behaviour in high-dimensional 
dynamical systems has gathered much attention []!]-§]. An interesting limit case of macro- 
scopic coherence is the global synchronization of all the parts of the system: this behaviour 
has been found in many different situations in physics, chemistry and biology. Relevant ex- 
amples of synchronization are the flickering of swarms of fireflies, fish locomotion, Josephson 
arrays and certain multi-mode laser systems |J. 

Beyond synchronization there exist much more interesting and intriguing phenomena 
involving non trivial coherence. Let us mention the relevant works of Kaneko showing 
the emergence of subtle coherence among different elements of globally coupled maps (GCM), 
in the case of a logistic local map. 

A GCM system (of the type considered in this paper) is a collection of elements x n 
(defined on the site n) evolving according to the equation: 

e N 

x n (t + 1) = (1 - e)f a (x n (t)) + - £ f a (Xi(t)), (1) 

iV i=i 

where N is the total number of elements, t is the discrete time and f a (x) is a nonlinear 
function depending on a parameter a. 

When the mapping function is f a (x) = ax(l — x) (i.e. a logistic map), for certain values 
of e and a one can observe rather remarkable behaviours, as anomalous fluctuations, glassy 
behaviour and broken ergodicity ||. For instance, Kaneko studied a parameters region 
where the number of positive Lyapunov exponents is proportional to iV; in this region the 
local elements x n take different values, almost randomly, so one might expect that the 
temporal fluctuations of the mean-field 

1 N 

fc(*) = ^E /.&(*)) (2) 

iV i=l 

obey a Gaussian statistics, with a mean square deviation a 2 N = (h 2 ) — (h) 2 proportional to 
N^ 1 . On the contrary one observes a 2 N — > const, as iV — ► oo. Only with the inclusion of 
noise one obtains a 2 N ~ iV _/3 , with (3 < 1, if the noise is not too strong. This is due to a 
coherence (collective behaviour) among the elements x n , that can be associated 0] with the 
existence of the periodic windows in the chaotic region of the logistic map. 
In some regions of the parameters space one observes clustering i.e., the elements x n 
split into clusters, such that the elements belonging to the same cluster k take an identical 
value Xk{t). The number C of the possible clusters can be very large with iV (InC ~ VN) 
and the dynamics can show glass-like behaviour and failure of ergodicity ||. 

Symplectic GCM, where the coupling is not of the mean-field type of eq. ([!]), can show a 
more intriguing behaviour for the clustered motion, since, in this case, the synchronization 
between the different parts is not perfect and the clustering has a finite lifetime ||. 

Even coupled map lattices (CML) with local coupling, on <i-dimensional lattices, can 
show a macroscopic behaviour that is very different from the microscopic one. For example, 
Chate and Manneville |J showed that in CML on hyper-cubic rf-dimensional lattices some 
macroscopic variables, as the center of mass: 
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1 N 

™(*) = T7 !>(*)> (3) 

i=l 

can have periodic behaviour, or quasi-periodic for d > 5, even if the system is chaotic, i.e. 
the first Lyapunov exponent is positive. 

In our opinion, one of the more interesting facts among the many different aspects of 
the high-dimensional dynamical systems is the macroscopic chaos. With this locution we 
indicate the irregular behaviour of some global variables (e.g. the center of mass) that evolve 
with typical times much longer than the characteristic time of the full dynamics, that we 
call the microscopic dynamics; the order of magnitude of the latter time may be estimated 
by the inverse of the first Lyapunov exponent Ai of the microscopic dynamics 

The macroscopic chaos for the high-dimensional maps, in some sense, is the analogous 
of the hydrodynamical chaos for the molecular motion. Let us briefly clarify this point and 
stress the open problems. For the Hamiltonian dynamical system describing the evolution of 
the molecules of a fluid the order of magnitude of the microscopic characteristic time is the 
mean collision time. On the other hand, if one studies the hydrodynamic equations of the 
considered fluid - those describing the motion on space and time scales that are very large 
with respect to the molecular scales - the involved characteristic times are much larger. So 
one expects that for the full Hamiltonian dynamics Ax ~ 10 11 s _1 , while the first Lyapunov 
exponent for the hydrodynamic equations, in a chaotic regime, is of order (10 _1 — 1) s _1 , 
the exact value depending on the geometry of the system, the external forcing, the Reynolds 
number and so on. 

This big difference comes from the coarse graining process in the building of the hy- 
drodynamic equations, that is from the very different scales one is observing in the two 
cases. 

If one is able to write down the equations at the hydrodynamical level (i.e. for the 
slow variables) then it is not difficult to characterize the macroscopic behaviour, by means 
of the standard techniques of dynamical system theory. However in this case, one subtle 
point, on the mathematical side, is the problem of the norm to be used in the definition 
of the Lyapunov exponents. For finite dimensional systems all the norms are equivalent 
|TjJ and there is no ambiguity. This is not true for infinite dimensional systems (as the 
hydrodynamical partial differential equations are), so that the Lyapunov exponents can 
depend on the used norm. 

In generic CML, as far as we know, there are no general systematic methods to build 
up the macroscopic equations; therefore it is necessary to study the macroscopic behaviour 
of the system in terms of the x n variables, i.e. it is necessary to remain at the microscopic 
level of description. 

The present paper is organized as follows. In the next section we report some phe- 
nomenology of macroscopic chaos in GCM. Section [III] is devoted to a discussion on the 
recently proposed self-consistent Perron-Frobenius method for the one-body description of 
these systems at a macroscopic level. In section [TV|, that is the nucleus of the paper, we show 
that the recently introduced finite size Lyapunov exponent gives a suitable quantitative char- 
acterization of the macroscopic chaos; that is to say: we are able to identify a macroscopic 
characteristic time by means of a Lyapunov exponent that is effective on the macroscopic 
scales and which, in general, is much smaller than the usually defined one. The last section 
is devoted to discussion and conclusions. In the Appendix we discuss some subtle points 
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related to the problem of the choice of the norm in the definition of the Lyapunov exponent 
for the one-body description. 

II. SOME PHENOMENOLOGICAL RESULTS AND FIRST NAIVE 

APPROACHES 

For sake of self-consistency, we report some numerical evidence (partly already known 
from previous works of the existence of macroscopic chaos in GCM, eq. ([I]). For a 

detailed description of the many different macroscopic behaviours see the works of Kaneko 
and also Pikovsky and Kurt lis 0. 

We studied mainly systems (|TJ) of globally coupled tent maps, i.e.: f a (x) = a(l/2 — 
\x — 1/2|), and all the figures of this paper refer to these systems. In addition, to test the 
validity of our method, we also considered a slightly different version of eq. (|l|), the so-called 
heterogeneous GCM J/J, in which the interacting elements are not all identical: 

e N 

x n (t + 1) = (1 - e)f an (x n (t)) + -£ /«,(*,(«)), (4) 

i=l 

the control parameters a* are quenched variables that take values according to some rule to 
be specified. Here we chose f ai (x) = 1 — a«x 2 (the logistic map), with a, assuming equally 
spaced values in an interval between a min and a max . 

At varying e and a in eq. (p]) or a min and a max , in eq. (f|), and by looking at the 
evolution of a macroscopic variable as, for instance, m(t), we observed different macroscopic 
behaviours, which can be grouped in the three following classes. 

i) STANDARD CHAOS: the microscopic Lyapunov exponent is positive and the law of 

large numbers holds. In this case, if one plots m(t) versus m(t — 1), one observes a 
spot with size of order 1/y/N. 

ii) MACROSCOPIC PERIODICITY: the microscopic Lyapunov exponent is positive and 

the plot of m{t) versus m(t — 1) gives rise to k spots with size of order 1/y/N, for a 
macroscopic motion of period k. There are also less trivial periodic or quasi-periodic 
behaviours (e.g. in a case is shown in which the mean-field dynamics evolves on a 
2-dimensional torus). 

iii) MACROSCOPIC CHAOS: the microscopic Lyapunov exponent is positive and the plot 

of m(t) versus m(t — 1) shows an apparently well structurated function, see fig. 1, that 
suggests chaotic motion for m(t). 

The first hope is that, at a macroscopic level, the system can be described by few global 
variables, that are able to catch the large scale properties, i.e. one hopes that the macroscopic 
motion is driven by a low dimensional dynamics. Indeed, fig. la (that is obtained for 
particular parameters values of the tent GCM), seems to suggest such a possibility: the 
first-return plot for the center of mass evolution appears to generate a one dimensional map. 
In such a case one can conjecture that 

m(t+l)=F(m(t)) + 7 lN (t), (5) 
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where t/n is a small (~ 1/vJV) incoherent fluctuation, vanishing in the limit iV — > oo. If 
the above equation holds - with the evolution function F(m) fitted from the data - one 
can define the Lyapunov exponent for dynamical system (0). This latter quantity gives 
information on the sensibility of the evolution of the macroscopic variable m on the initial 
conditions, so that one may call it macroscopic Lyapunov exponent A TOacro ; according to the 
standard definition, one has: 

^macro = / ln \F'(m)\P(m)dm, (6) 

where P(m) is the probability distribution of m, as obtained from numerical simulations of 
eq. (0). In the case of macroscopic chaos one expects that A macro > 0. On the contrary, 
using the recipe of eq. (|6|) one obtains that A macro is negative. This result is not a paradox: 
simply eq. (|) does not hold. If one looks carefully at the plot of m{t) versus m(t — 1), 
for large N, one realizes that very fine grained structures exist, see fig. lb. This suggests 
that the variable m(t) alone is not sufficient for a satisfying description of the macroscopic 
properties of the system. Therefore we need to introduce some other variables. It seems 
likely that the system evolution at a macroscopic scale is driven by a full set of equations of 
phenomenological type, such as, e.g.: 



m 
a 



(t + l) = F 1 {m{t),a{t),...)+r ] ^ 

(t+1) = F 2 (m(t),a(t),...)+ V ^ (7) 



where a 2 = {J2x 2 )/N — (J2xi/N) 2 , and T)® are small fluctuations vanishing as N — > oo. 
Unfortunately, as far as we know, there is not a systematic approach (even of numerical 
type) for the building of eq.s (0), i.e. for the determination of the functions F 1; F 2 , .... 

Another naive approach, based on the possible time scales separation between the mi- 
croscopic and the macroscopic dynamics, is the analysis of the microscopic and macroscopic 
time correlation functions, defined as follows: 

r {i) / x {Xjjt + T)Xi{t)) - (Xj) 2 

^microV) / 2\ _ / X -\2 ^ ' 



_ (m(t + T)m(t)) - (m) 2 

^macroXj ) — / 2\ / \2 ' V / 

{m z ) — (m) 2 

where ((...)) means time average along the trajectory. 

It is to be noted, at this point, that the system we are studying shows a non trivial 
macroscopic behaviour in a parameter region where it is also in a broken ergodicity phase. 
In fig. 2 we show the probability distributions, P n (x), of two different x n for the system ([[]) 
in a situation of macroscopic chaos. It is well evident that the knowledge of the probability 
distribution of a generic x n does not permit a good description of the system (this is what 
we mean, in such a setting, by saying that the system is not ergodic). It is important to 
note that the two elements have been chosen in the two different bands in which the {x n } 
are grouped, see [EJ. Moreover elements belonging to the same band seems to have the 
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same distributions. Let us stress that this behaviour is not related to the particular band 
structure of the tent GCM: we found non trivial macroscopic motion in the presence of 
broken ergodicity also in heterogeneous GCM ([|), where the band structure is absent, but 
the microscopic building elements are not identical. 

Fig. 3 shows the above defined correlation functions in the regime of macroscopic chaos. 
Since one has a failure of the ergodicity, different sites have different autocorrelations, that 
is Cmicro{ T ) depends on the band to which the element X{ belongs. Rather unexpected is 
the fact that the characteristic decay times of C,^ cro (r), T l m , and of C macro (r), tm, defined 
by the equations 

CLcroir) ~ e-^, C macro (r) ~ , (10) 

are not very different; for instance, in the case of fig. 3 one has: T % m « tm ~ 115. A non 
trivial point is that these typical decay times are much larger than the characteristic time 
given by the Lyapunov exponent, i.e. 1/Ai (for the system of fig. 3 one has 1/Ai ~ 4.). 
This fact can be related to the strong intermittency of the system. We have to stress that 
this behaviour of the correlation times has been observed also in the heterogeneous logistic 
GCM (§. 

In the next sections we discuss some less naive approaches to the characterization of the 
macroscopic chaos. 



III. ONE-BODY DESCRIPTION IN TERMS OF A SUITABLE 
PERRON-FROBENIUS OPERATOR 

Let us now discuss how one can try to characterize the coherent dynamics of GCM in 
terms of a suitable Perron- Frobenius (PF) operator |2|,|6|||. At any time t one can wonder 
about the snapshot density pt(y), defined by the probability density that x n (t) assumes the 
value y: 

1 N 

= ■*,£*(!/-**(*))■ (ii) 

iv i=l 

Using the fact that, for a mean field coupling as in eq. ([!]), each x n interacts with all the 
other {xi} in the same way, the time evolution of pt{y) can be written in terms of pt-i(y), 
by means of a self-consistent approach: 

. w . £ a Vt _ lW ._l_ ? g||, (12) 

where the sum is on the preimages given by 

y, = f^( y - tU f^- l(z)dz ), (13) 

f^ 1 being the inverse function of f a , for instance, if f a (y) = 1 — ay 2 then f^ 1 {z) = 
±y/(l - z)/a. 
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In the parameter region of the standard chaos (see sect. 0) where m(t) = const. + 
0(1/ VN), one has that, after long times, pt(y) converges toward a stationary distribution 
Peq(y), i-e.: 

Pt(y) ~> Pe q (y) for t^oo. (14) 

On the other hand, for a macroscopic periodic behaviour, the long time properties of the 
density ptiy) are described by means of a finite collection of functions of y\ e.g. when m{t) 
asymptotically oscillates among two different values: 

m(2t)=m e m(2t + 1) = m OJ (15) 

after a transient time one observes 



P2t{y) = p e (y) P2t+i(y) = Po(y)- (16) 

Finally, in the case of macroscopic chaos Pt{y) does not converge toward any finite set of 
function of y, but evolves following a nontrivial time behaviour. 

The above approach for pt{y) has been introduced by Perez and Cerdeira |J and K. 
Kaneko ||. From a conceptual point of view eq.s ([12] - [13]) are close to the Boltzmann 
equation for the gases ||11|| . Indeed, as one can check by a direct inspection, the PF operator 
£( N ) that rules the evolution of the probability density pi N \yi, y 2 , • • • , 2//v) i n the whole 
phase space: 



Pi 



(N) 



(xi,x 2 , • • .,x N ) = £ {N) pi-l{xi,x 2 , . . .,x N ), (17) 



is a linear operator. On the other hand the PF operator £W for the evolution of Pt(y), given 
by eq.s ([12] - [Of), is a nonlinear one ||, since it depends explicitly on In the kinetic 



theory of gases one has a similar situation: the Liouville equation, for the evolution of the 
probability density in the whole phase space, is linear, while the Boltzmann equation for the 
one-body statistics is nonlinear. 



Since the evolution equations (|T^ - 13) for pt(y) are a dynamical system (though of 
infinite dimension), one can formally introduce the notion of first Lyapunov exponent 0. 
By so doing, as noted in the introduction, one is faced with the mathematical problem of 
the choice of the norm: for sake of simplicity we adopt the usual norm in the L 2 space. One 
could think that the problem of the norm is relevant just for strictly mathematical reasons, 
while from a physical point of view this problem practically does not exist. In the Appendix 
we discuss how even for very simple shapes of p(x) one obtains different Lyapunov exponents 
using different norms. 

However to practically obtain the first Lyapunov exponent, firstly one applies an in- 
finitesimal perturbation 5po(y) on pt{y) (in order to satisfy the normalization condition one 
has to impose the constraint / Sp(y)dy = 0) then one observes the evolution Sp n (y) of the 
perturbation. In this way, following the standard methods of Benettin et al. [|l|], it is 
possible to estimate the whole set of the Lyapunov exponents {Aj}; in particular one has 

A, = 1 li m hmi£\n%±« (18) 
2T^oo 5 p -oT^ f\6p n (y)\ 2 dy K } 



7 



It is easy to understand that Ai can give interesting "global information" which is com- 
pletely different from the one given by the microscopic Lyapunov exponent. In the case 
of standard chaos, or macroscopic periodic behaviour, Ai is negative, or zero, while for 
macroscopic chaos Ai > 0. Let us stress that Ai can be negative even if the first Lyapunov 
exponent for the dynamical equations (|l|) is positive. Some detailed numerical investigation 
of Kaneko on GCM show that Ai is positive in the cases of macroscopic chaos f2|. 

Equations ([12] - ^) hold only in the case of a mean field interaction for homogeneous 
GCM (eq. |l|) and other particular systems ||, but they are not valid even for the heteroge- 
neous GCM (eq. [|). Nevertheless, the computation of Ai according to the definition (|TS|), 
remains possible using a numerical approach - at least in principle, even if it may result 
in a very huge and heavy task - also for generic CML (i.e., in the case of maps with local 
coupling). This is accomplished by following the evolution of pt(y), as given by the time de- 
velopment of the histogram of the instantaneous microscopic configuration (x%, X2, ■ ■ ■ , xn). 

Besides the previous relevant properties of Ai, we have to underline a delicate point that 
rises for its numerical computation, when it is not possible to write down the appropriate 
eq. (|T2|). In this case, as we already pointed, one is forced to study the evolution of \\5p t (y) \ \ 
making histograms from two instantaneous configurations {x{\ and {a^}, obtained by the 
introduction of a small perturbation on {xj} at t = 0. To show where is the problem, let 
us consider the time evolution of Sm(t) in the case of standard chaos. One can write 5m in 
two different ways: 

1 N 

Sm(t)^-J25x n (t) (19) 

iV n=l 



5m(t) = J y5 Pt (y)dy. (20) 

If we are in the parameters region where the Lyapunov exponent Ai is positive, but Ai is 
negative, we have both |5x n (i)| ~ |5a; n (0)| exp(Ait) and \\5p t \\ ~ | |<5po| I ex P( — |Ai|t). So that 
apparently we fall into a contradiction, since from flT9| ) we obtain 

\Sm{t)\ ~e Al ', (21) 

while from (p0|) we have 

\5m(t)\ ~ e- |Xl1 *. (22) 

To eliminate the contradiction one needs to take into proper account the scales of the 
perturbations on which these two equations hold true; eq. (|2lT) holds at small scales: \Sm\ <C 
0(l/y/N), while eq. (|2|) holds for \Sm\ > 0(1/ VN). Basically if one starts with a 
perturbation 5m(0) much smaller than 0(1/ VN) one has an exponential growth, but only 
up to a magnitude 0(l/y/N), the maximum allowed by the laws of large numbers. In an 
analogous manner starting from a perturbation 5m(0) much larger than 0(1/ ^/N) one has 
an exponential decrease up to 0(l/y/N), see fig. 4. 

From the previous example it is clear that if one studies numerically the time evolution of 
\\5pt\ \ for large but finite N, at small times (or, more precisely, small ||5p||) one can observe 
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just the "microscopic" Lyapunov exponent. So one has to introduce some coarse-graining - 
i.e. to consider non infinitesimal 8po - to obtain a macroscopic level description. However, 
it must be noted that there is a large arbitrariness in choosing a non infinitesimal Spo. 

In this section we have discussed the computation of the macroscopic Lyapunov exponent 
by means of the Perron- Frobenius method, both in its original formulation eqs. ([3~2| [3~3|) and 
in its numerical version, to be applied when the P-F equation cannot be explicitly written. 
Neither approach is easy to handle because of the arbitrariness in the choice of the norm 
and of the non infinitesimal Spo. 

In order to avoid these difficulties (numerical and conceptual) concerning the one-body 
description, in the next section we introduce a different approach. 



IV. CHARACTERIZATION OF THE MACROSCOPIC CHAOS VIA THE 
PREDICTABILITY IN THE LARGE 



The introduction of efficient methods (Benettin et al. |T^]) for the computation of the 
Lyapunov exponents has been an important milestone in the understanding of chaotic sys- 
tems. However, in spite of their deep conceptual relevance, the Lyapunov exponents are not 
able to give a complete characterization of realistic dynamical systems, such as systems with 
intermittent behaviour or with many different characteristic times. For example, in order 
to describe intermittency effects one needs to introduce generalized Lyapunov exponents 



and the statistics of the fluctuations of the effective Lyapunov exponent ||13|-|15|1. On the 
other hand in systems with many different characteristic times (e.g. the fully developed 
turbulence) the Lyapunov exponent cannot be used to study, even at very rough level, the 
predictability problem for small (but not infinitesimal) perturbations, that cannot be de- 
scribed in terms of the tangent vector. In order to overcome these difficulties recently a new 
indicator has been introduced: the finite size Lyapunov exponent (FSLE), see Ref. JTB| for 
details. 

The definition of FSLE may be given in terms of the predictability time T r (S), that 
is the time a perturbation of initial size 5 takes to grow by a factor r (> 1) during the 
system evolution. The perturbation of size 5 is supposed to be already aligned with the 
most unstable direction. The factor r should be taken not too large, in order to avoid the 
growth trough different scales. In many applications, r = 2, so that sometimes the T r (8) is 
also called the error doubling time. The Finite Size Lyapunov Exponent is defined from a 
suitable average of the predictability time according to 



>(i) iwr-(w hr (23) 

where (• • denotes average with the natural measure along the trajectory and (• • ■) is the 



average over many realizations, for more details see ref. [H]]. The above definition holds 
for continuous time system. In the case of discrete time systems (i.e. maps) T r (S) is the 
minimum time such that the distance A r between the two realizations is greater or equal to 



r5. Therefore, instead of eq. fl23|) we have \T%: 



m= <mWT»- (24) 
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Let us note that by the definition ( p3| ) A (5) is not appropriate in discerning cases with 
A = and A < 0, due to the positiveness of the predictability time. 

For chaotic systems the limit of infinitesimal perturbations 5 — > in ( p3|) gives the 
maximum Lyapunov exponent Ai, i.e. A (5) displays a plateau at the value Ai for sufficiently 
small 5. 

In many realistic situations, the error growth for infinitesimal perturbations is dominated 
by the faster scales, that are typically the smaller ones (as in the classic example of three 
dimensional turbulent flows). When the size 5 of the perturbation cannot be considered 
any longer infinitesimal, the behaviour of A (5) is governed by the nonlinear evolution of 
the perturbation, and, in general, X(5) < X\. The decrease of X(5) does follow a system- 
dependent law. In some cases, A (5) can be predicted by dimensional arguments, e.g. in the 
fully developed turbulence one has the universal law A (5) ~ 5~ 2 , in the inertial range ||16| , |T7 



Therefore, the behaviour of A as a function of 5 gives important information on the 
characteristic times (and scales governing the system and is a powerful tool in studying 
dynamical systems which involve many characteristic scales in space and time. 

The FSLE has been successfully used even in time series analysis [flD ], A relevant practical 
fact is that, at variance with the usual Lyapunov exponent which needs high embedding 
dimensions and extremely huge statistic, the FSLE at relatively large 5 can be extracted even 
using moderately small embedding dimension. For example, in systems in which some slow 
variables can be separated from the fast variables, the FSLE related to the slow dynamics 
can be obtained even with moderate statistics and unresolved small scales (i.e. the fast 
ones) . 

To practically compute the FSLE, firstly one has to define a series of thresholds 5 n = r n 5o 
in = 1, . . . , M), and to measure the times T r (5 n ) that a perturbation of size 5 n takes to grow 
up to 5 n+ i. The times T r (S n ) are obtained by following the evolution of a perturbation from 
its initial size 5 m i n < 5o up to the largest threshold 5m- This can be done by integrating 
two trajectories of the system starting at an initial distance 5 m i n . The FSLE, \(5 n ), is 
then computed by averaging the predictability times over several realizations, see equation 
(^U). In general, one must choose 5 m i n <C 5q, in order to allow the direction of the initial 
perturbation to align with the most unstable direction in the phase-space. Moreover, one 
must pay attention to keep 5 m < 5 satU rationi so that all the thresholds can be attained. 

Our hope is that, as in the usual hydrodynamical description of thermodynamic systems, 
the variables of systems with macroscopic chaos can be separated in slow variables (given 
by m(t), a(t) and so on) and fast variables. 

We have computed the FSLE fl23|) looking only at m(t) both for the homogeneous GCM 
(0), using as local map the tent map, and for the heterogeneous GCM (B), for different 
values of the parameters. Practically the numerical experiments has been performed as 
follows: we have averaged, over many realizations, the predictability times T r (5 n ) (defined 
above) related to |5m(t)| which is initialized at the value 5m{t) = 5 min . This has been done 
considering, for each realization, the perturbed system that is obtained by shifting all the 
elements of the unperturbed one by the quantity 5 min (i.e. £j(0) = x«(0) + 5 min ). 

As can be seen from figure 5a, which shows A (5) versus 5 in the case with macroscopic 
chaos, two plateaux exist: at small values of 5 (5 < 5i) one has, as expected from general 
considerations, A (5) = Ai; for 5 > 52 one has another plateau which can fairly be called 
"macroscopic" Lyapunov exponent, A (5) = X macro . Moreover, 5i and 5 2 decrease at increas- 
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ing N. In fact, from figure 5b, where we report A (5) versus SyN for different values of 
N, one can see that both 5\ and 5 2 scale as ~ \/v~N\ by rescaling 5 with vN, one ob- 
serves that SivN (which indicates the end of the microscopic plateau) and 52\f~N (which 
marks the end of the saturation due to the microscopic dynamics and the beginning of the 
macroscopic plateau) coincide for all values of N. Incidentally, one can also observe that 
the macroscopic plateau, being almost non-existent for N = 10 4 , becomes more and more 
resolved and extended on large values of 5yN at increasing N up to N = 10 7 . From these 
features we can argue that the macroscopic motion is well defined in the limit N — > oo and 
therefore we can conjecture that in this limit the microscopic signature in the evolution of 
Sm(t) completely disappears in favor of the macroscopic behaviour. 

It is important to note that in order to obtain the results shown in fig. 5 and fig. 7 some 
precautions must be taken in initializing the system. In fact, because of the complex struc- 
ture of the phase space of this system, one has that different initial conditions can led to 
different attractors ||. Therefore, choosing random initial conditions while increasing N, 
one likely obtains different values for the macroscopic and/or microscopic Lyapunov expo- 
nents. In order to maintain the system on the same attractor, we kept fixed the percentage 
of elements per band (see Kaneko @), at varying the number of elements N. 

In the case of standard chaos (A macro < 0) one has only the microscopic plateau and then 
a fast decreasing of \(S), as can be seen in fig. 6. Also in the cases of non trivial macroscopic 
behaviour, either periodic or quasi-periodic (i.e. A macro = 0), we see only the microscopic 
plateau followed by a decreasing up to the saturation. This is a consequence of the fact that 
X(S) is not able to distinguish between A < and A = 0, by definition, see eq. (]23). 

One can wonder if it is possible to compute the macroscopic Lyapunov exponent directly 
looking at (In \5m\) versus t. The macroscopic regime is expected to appear as a change of 
slope when the size of the perturbation is 0(1/ yN) (as discussed in sect. |TTTj in the case 
of standard chaos). In fig. 7 (In \ Sm\) is shown as a function of t for the same systems of 
fig. 5 (with identical initial conditions). One can see that for small values of 5m(t) one has 
a well defined exponential growth according to the microscopic Lyapunov exponent. This 
regime stops at 8m(t) ~ 0(1/ ^/N) i.e. when the microscopic error growth saturates. At 
larger values of Sm(t) one observes a second (macroscopic) slope. One has a clear evidence 
also from fig. 7 that the macroscopic motion sets in at a scale that is smaller when N is 
larger, as already established by the analysis of fig. 5. In the case of fig. 7 one has also a fair 
quantitative agreement with the results obtained with the FSLE method (fig. 5). However 
in general, because of the intermittency, the macroscopic regime can start at different times 
and with different rates, therefore the naive method in which the averages are performed at 
fixed delay times can give spurious laws; see ref TOJIH | for a discussion of these effects. 

Of course a priori there is no deep reason to believe that the above analysis performed 
only on one variable, m(t), is enough accurate, therefore we performed a similar study 
looking at two macroscopic variables, m(t) and o~(t). Since there is not a natural way to 
introduce a perturbation on m(t) and cr(t), we are forced to follow the method of Wolf et al. 
PHI for the computation of the Lyapunov exponent from data, but with some modifications. 
We considered the phase space generated by the variables m(t), a(t) and we looked at the 
time separation between analogues whose initial distance 8q is not very small (i.e. in the 
region of the "macroscopic" plateau in fig. 5). Following the evolution of the full system we 
compute the quantity: 
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where A(t) = ^5m 2 (t) + Sa 2 (t), tj are the "return" times such that A(t,) e [5o,2 5 ] an d A 
is the number of analogues. Fig. 8 shows T(r) versus r. 

It is important to stress the relevance of the value of Sq, which must be not too small. 
Using a 5 out of the "macroscopic" (i.e. 5 < 5 2 ) plateau one has spurious behaviour due 
to intermittency (see ref. flSH ). 

It is remarkable the agreement between the values for the macroscopic Lyapunov expo- 
nent obtained with the two methods. 



V. CONCLUSION AND DISCUSSION 

In this paper we proposed a quantitative characterization of the coherent evolution in 
globally coupled maps that show macroscopic chaos. From a conceptual point of view, the 
irregular behaviours of global observables - which typically evolve on time scales much longer 
than those related to the evolution of a single element x n (t) - are for the GCM the analogous 
of the chaos for hydrodynamical equations in the description of the molecular motion. 

In generic systems there are not a systematic method to obtain the equations that de- 
scribe the evolution of the degrees of freedom accounting for the large scale behaviour. 
Therefore a straightforward application of the usual dynamical systems methods is not able 
to give suitable characterizations. In particular, the usual Lyapunov exponent, which gives 
a microscopic description, is not able to capture the essential features of this macroscopic 
phenomenon. 

It is worth to note that in the cases where it is possible to derive hydrodynamical-like 
equations, such as the Perron-Frobenius equation, one obtains a description of the system 
at a macroscopic level, see for a detailed discussion @||. We have shown that, in the cases 
in which this approach is missing, the recently introduced Finite Size Lyapunov Exponent 
T6| is able to characterize, in a consistent way, the macroscopic behaviour of the system 



and gives a measure of the chaotic motion at a coarse-grained level. We have the following 
scale-dependent scenario for the mechanism acting on the evolution of the perturbation of 
global observables (i.e. the mean field): 



at small values of the perturbation (<C 1/yN) we recover the microscopic Lyapunov 
exponent, simply because at these scales the evolution of the perturbation is driven 
by the microscopic dynamics; 



• at larger values (^> 1/yN) we obtain a macroscopic Lyapunov exponent, which can 
be much smaller than the first Lyapunov exponent of the microscopic dynamics. 

The above scenario is confirmed using (with some important changes) the method of Wolf 
et al. [2(| for the computation of the Lyapunov exponent from numerical data. Basically, in 



presence of macroscopic chaos one observes that: the variables of the system are separated 
in slow variables (e.g. m(t),a(t)) and fast variables (see fig.s 5 and 7). For a discussion of 
the FSLE in systems with slow and fast dynamics see Ref. [plf] . 
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Let us stress that also with a naive numerical treatment of the one-body description one 
can observe a rather clear evidence for the two regimes (the microscopic and the macroscopic 
ones) characterized by two very different "Lyapunov exponent". But, as discussed in the 
Appendix, the problem of the choice of the norm (that is present also working with the 
histograms) gives some unavoidable troubles which do not allow for a detailed comparison 
with the other methods. 

Finally, it is interesting to mention the case of "single-band" chaos (see ref. ||), where 
one has irregular behaviour with some macroscopic coherence, but there is not a separation 
between the microscopic and the macroscopic scales. The microscopic Lyapunov exponent 
has the same value of the macroscopic one (A macro = Ai). This means that in this case - as 
fig. 9 shows - A (5) has a single plateau, extending over a range of scales larger than that of 
the standard chaos (fig. 6), so as to include the macroscopic scales. Let us remark that this 



case of "single-band" chaos is reminiscent of the two-dimensional turbulence [H]] , where the 
A (5) does not depend on 5 in a very large range of values; that is, a broad range of scales 
have equal characteristic times. 



After the submission of this paper, the authors are informed of the recent preprint by T. 
Shibata and K. Kaneko "Collective Chaos" ( |chao-dyn/9805009|) , where a related study is 



presented. In particular, T. Shibata and K. Kaneko use the FSLE to study the macroscopic 
properties of a heterogeneous GCM. Their results are in agreement with ours. 
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APPENDIX: 

One could naively think that the behaviour of |#m(i)| and \ \Spt\ \ give the same informa- 
tion. In this appendix we show that the above idea does not work because of the unavoidable 
troubles related to the choice of a norm for Sp t . 

In order to show this we discuss two explicit cases which clarify this point. As first 
example we consider p(x) uniform in the interval [a, b] and the perturbed distribution p (x) 
uniform in the interval [a + 5m, b + 5m], i.e. with a shifted average. Let us now compute 
two different norms of 5p = p — p: 

/ \6p{x)\dx = -^-\6m\ (Al) 
J b — a 




\\5p\\ 2 = J / 8p\x)dx = J— -\5m\ . (A2) 
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Assuming that 8m is related to a GCM system, it is clear that the computation of the 
Lyapunov exponents according to the norms defined above gives, in this case (uniform 
distribution), two Lyapunov exponent different by a factor two. Nevertheless with different 
shape of the distribution we have different results. To show this last issue we consider now 
a p(x) with a symmetric triangular shape in the interval [a, b] and zero for x < a or x > b. 
The perturbed p has the same functional form with the average shifted of 8m. The same 
computations as before lead to: 



From the above examples it is clear that the relation between \8m\ and \\8p\ \ depends 
dramatically on the explicit norm used and on the shape of the p. In addition, in the real 
time evolution of 8p t one has not only shifting of the average but also variations of a 2 and 
other parameters of the distribution that complicate the interpretation of \\8p\ \ in terms of 
\8m\. Therefore there is no reason to expect a quantitative agreement between numerical 
study (using histograms) of the one-body description and \8m(t)\. 




(A3) 




(A4) 
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FIGURE CAPTIONS 



FIGURE 1: (a) Plot of m(t) versus m(t — 1) for the system (|]) with a = 1.7 e = 0.3 and 
iV = 10 6 ; fig (lb) is an enlargement of (a). 

FIGURE 2: Probability distribution function P n (x) for two different x n for the system eq. 
(|TJ) with a = 1.8, e = 0.3 and iV = 10 5 . The distribution are computed following x n (t) 
for 5 • 10 5 temporal steps. 

FIGURE 3: Absolute value of the correlation function C macro {r) (solid line) compared 
with the absolute value of the correlation functions for two different Xj(t), C l micro {r), 
(dashed and dot-dashed lines). The parameters are chosen as in fig. 2. 

FIGURE 4: < In |<Jm| > versus t for the system (0) with a = 1.8, e = 0.08 and N = 10 4 . 
The squares are for an initial perturbation 5m(0) = 10~ 6 <C 1/vA, the circles for 
5m(0) = 0.15 ^> 1/y/N, and the dotted line has slope X m i cro . The average is over 10 3 
realizations. 

FIGURE 5a: X(5) versus 5 for the system flip with a = 1.7, e = 0.3 for N = 10 4 (x), 
N = 10 5 (□), N = 10 6 (0) and N = 10 7 (A) . The first plateau corresponds to the 
microscopic Lyapunov exponent X m i cro ~ 0.17 and the second one to the macroscopic 
Lyapunov exponent A macro ~ 0.007. The average is over 2 ■ 10 3 realizations for iV = 
10 4 , 10 5 , 10 6 and 250 realizations for iV = 10 7 . The system has been initialized in 
order to leave the attractor unchanged (as explained in the text). 

FIGURE 5b: The same as figure 5a rescaling the S— axis with v^A- 

FIGURE 6: X(S) versus 5 for the system (0) with a = 1.8, e = 0.08 for N = 10 4 (□), 
iV = 10 5 (©) and iV = 10 6 (A). The plateau corresponds to the microscopic Lyapunov 
exponent A m j cro ~ 0.5. The average is over 2 • 10 3 realizations. 

FIGURE 7: < In \Sm\ > versus t for the system with a = 1.7, e = 0.3 and N = 10 4 
(upper curve), iV = 10 5 (middle curve), = 10 6 (lower curve) The straight lines 
correspond to the two exponential growth the microscopic one (A micro ~ 0.17) and the 
macroscopic (A macro ~ 0.01). The average is over 500 realizations. 

FIGURE 8: T(t) versus r for the system (0) with a = 1.7, e = 0.3 and A = 250, the 
dashed line corresponds to A macro ~ 0.006. The distance has been plotted every two 
time steps to eliminate oscillations given by the swiching between bands. 

FIGURE 9: X(8) versus 5 for the system (g) with a = 1.5, e = 0.3 for N = 10 4 (□), 
iV = 10 5 (©). The plateau corresponds to a Lyapunov exponent X m i cro ~ 0.28. The 
average is over 2 • 10 3 realizations. 
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